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Abstract 

We propose an infrared cut-off for the holographic density which incorporates 
among others a constant term, that produces the effect of the cosmological 
constant, improving the results of previously considered holographic models 
based on local quantities. The inclusion of constant term is interpreted as a 
natural first approximation for the infrared cutoff which is associated with the 
vacuum energy, and the additional terms guarantee an appropriate evolutionary 
scenario that fits the astrophysical observations. Cosmological constraints on 
the model have been studied using the observational data. 

PACS: 98.80.-k, 95.36,+x 

1 Introduction 

The astrophysical data from distant la supernovae observations [T], [2J, [3J, [I], [5], IS], 
[7], [8], cosmic microwave background anisotropy [9], pU], [H], [12], and large scale 
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galaxy surveys [13], [H], [15], [16], [IT], [18], [19], all indicate that the current Uni- 
verse is not only expanding, it is accelerating due to some kind of negative-pressure 
form of matter known as dark energy ([2D], |21j.|22].|23]. |24j . |25j). The combined 
analysis of cosmological observations also suggests that the universe is spatially flat, 
and consists of about ~ 1/3 of dark matter (the known baryonic and nonbaryonic 
dark matter), distributed in clustered structures (galaxies, clusters of galaxies, etc.) 
and ~ 2/3 of homogeneously distributed (unclustered) dark energy with negative 
pressure. Despite the high percentage of the dark energy component, its nature as 
well as its cosmological origin remain unknown at present and a wide variety of mod- 
els have been proposed to explain the nature of the dark energy and the accelerated 
expansion (see [2D]-[25] for review). Among the different models of dark energy, the 
holographic dark energy approach is quite interesting as it incorporates some concepts 
of the quantum gravity known as the holographic principle ([261 123 123 129| 130] ) .which 
first appeared in the context of black holes [21] and later extended by Susskind [3D] 
to string theory. According to the holographic principle, the entropy of a system 
scales not with its volume, but with its surface area. In the cosmological context, the 
holographic principle will set an upper bound on the entropy of the universe [31] . In 
the work [2S], it was suggested that in quantum field theory a short distance cut-off A 
is related to a long distance cut-off (infra-red cut-off L) due to the limit set by black 
hole formation, namely, if is the quantum zero-point energy density caused by a short 
distance (UV) cut-off, the total energy in a region of size L should not exceed the 
mass of a black hole of the same size, thus L 3 A 4 < LMp. Applied to the dark energy 
issue, if we take the whole universe into account, then the vacuum energy related to 
this holographic principle is viewed as dark energy, usually called holographic dark 
energy [2D] [22], [S3]- The largest L allowed is the one saturating this inequality so 
that the holographic dark energy density is defined by the equality pa = 3c 2 M 2 L~ 2 , 
where c 2 is a numerical constant and M~ 2 = 8ttG. 

As is well known, the Hubble horizon H~ x as the infrared cut-off, gives an equation 
of state parameter (EoS) equal to zero, behaving as pressureless matter which cannot 
give accelerated expansion, and the particle horizon gives an EoS parameter larger 
than —1/3, which is not enough to satisfy the current observational data. The in- 
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frared cut-off given by the future event horizon yields the desired result of accelerated 
expansion with an EoS parameter less than —1/3, despite the fact that it is based on 
non local quantities and has problems with the causality [33]. An holographic dark 
energy model which is based on local and non local quantities have been considered 
in [3l], [35]. Based mainly on dimensional arguments, in |36] . [37] an infrared cut-off 
for the holographic density was proposed, which describes the dark energy in good 
agreement with the astrophysical data, and may explain the cosmic coincidence. This 
model exhibits quintom nature (for some values of the parameters) without the need 
to introduce any exotic matter, and has also proven to be useful in the reconstruction 
of different scalar field models of dark energy which reproduce the late time cosmo- 
logical dynamics in a way consistent with the observations [38] [39] . 
Despite the available proposals, the theoretical root of the holographic dark energy is 
still unknown mostly because our ignorance about the microscopic nature of quantum 
gravity and the origin of dark energy. This translates into a lack of details about the 
holographic cut-off, and therefore in the present paper we assume that this cut-off 
should be encoded in a function that depends on natural cosmological quantities like 
the Hubble parameter, its time derivative and also should contain a constant term, 
which represents the first approximation to vacuum energy density (which could be 
interpreted, in the context of the holographic principle as giving quantum gravity na- 
ture to the cosmological constant). The constant term would act as the first natural 
approximation to the infrared cut-off for the holographic density, since this term may 
be identified with a constant vacuum energy density, that adjusts very well to the 
current observational data. A discussion about the cosmological constant and the 
holographic principle is given in [ID]. We propose an holographic density given by a 
general function of the form f(H,H), where / should be expanded in powers of his 
arguments, based on their smallness (compared to the squared Planck mass M£). The 
coefficients of the different terms of this expansion must be found by fitting the model 
with current observations. To obtain the cosmological constraints on the model, we 
will use the SNIa Constitution sample ([H]), complemented with the CMB (cosmic 
micro wave background) anisotropy and BAO (baryon acoustic oscillation) data. In 
section 2 we review the main aspects of the model, in section 3 we include the matter 
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content and fit the parameters using the joint Constitution, CMB and BAO data 
analysis. In section 4 we give a discussion. 



2 The Model 

The main inequality discussed above L 3 A 4 < LM 2 , leads to the restriction involving 
the UV and IR cut-offs A 4 < M 2 L~ 2 . Saturating this inequality we come to the 
concept of holographic dark energy p\. To define the infrared cut-off and hence the 
holographic dark energy density, we propose for L~ 2 a general function which should 
depend on the main cosmological parameters H and H, and is measured in units 
of (length)~ 2 , defining the holographic density as follows (from now on we will set 
8nG = M- 2 = 1) 

p A = 3f(H 2 ,H) (2.1) 

where H = a/a is the Hubble parameter, and the function / can be Taylor expanded 
in powers of his arguments provided H 2 and H are at the present very small quantities 
measured in Planck mass units, giving 

f(H 2 ,H) = X + aH 2 + f3H (2.2) 

where we neglected higher order powers of H 2 and H, and A, a and /3 are constants to 
be determined. This model generalizes the one presented in (|36].|37]). The novelty 
of the present proposal takes root in the presence of the constant term A which 
we consider as the first natural approximation to the infrared cut-off (L ~ A -1 / 2 ), 
and should be associated with the contribution of the constant energy density of 
the vacuum. As we shall see bellow, the cosmological dynamics depends on a and 
(3 only through the combination (a — l)//3 (i.e. cosmological observations fix only 
this combination) leaving one free parameter. To fix the free parameter we can 



use geometrical criteria in order to reduce the last two terms in (2.2) to a term 
proportional to Ricci scalar, but the dynamics is independent of the election of this 
free parameter. 



To analyze some cosmological consequences of the model (2.2), we first consider the 
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case of an universe dominated by dark energy. The Friedmann equation in absence 
of matter is 

H 2 = \p A = A + aH 2 + (3H (2.3) 



Setting x = In a, we can rewrite (2.3) as follows 



H i = X + aH 2 + ?- d ^- (2.4) 
2 ax 

Introducing the scaled Hubble expansion rate H = H/Hq, where Ho is the present 
value of the Hubble parameter (for x — 0), the above Friedman equation becomes 

H 2 = Q A + aH 2 + ^ d ^- (2.5) 
2 ax 

where Q A = X/Hq Solving the equation (2.5), we obtain 



H 2 = Q A + VLe- 2x{a ~ 1)/p = Vt A + n{l + z) 2{a - 1)/p (2.6) 

where Q A = Q\/ (1 — a), is an integration constant and the last equality is written 
in the redshift variable (by using e~ x = 1 + z). Note that even without the matter 
component, this solution is interesting enough as it contains the cosmological con- 
stant, and the second term may give matter-like behavior if the constants a and /3 
satisfy the restriction (a — l)//3 = 3/2, giving rise to the ACDM model. Evaluating 



the Eq. (2.6) at x = (the flatness condition), it follows 

n A + n = i. (2.7) 

Thus f^A represents the constant component of vacuum energy density and Q stands 
for the new form of "matter", both quantities coming from the holographic principle. 
Additional constraints may be used to find appropriate values for the parameters. 
These constraints are related with the present observational value of the dark energy 
EoS, which between the margins of errors can be set to w ~ —1.1, and the transition 
deceleration-acceleration, which occurs in the region zt ~ 0.5—0.7 (zt is the transition 
redshift). The EoS in terms of the redshift is given by 

1(1 + z)d(H 2 ) , . 

w(z) = -l + - v ~ ' \ J 2.8 
K J 3 H 2 dz K J 
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replacing H 2 from (2.6) it follows that 



w(z) 



1 - 



2(1 - a) (I - fl A )(l + zf^- 1 ^ 



(2.9) 



(2.10) 



3/3 [n A + (l-fi A )(l + z) 2 («-D// 3 ] 
from Eq. (2.6) the deceleration parameter q = —ad/a 2 , is given by 

n (l-a)(l-n A )(l + *) 2(a - 1)//? 

q[Z) ~ /3[fi A + (l-fi A )(l + z) 2 (- 1 )// 3 ] 

Note that the Hubble parameter ( 2.6[ ) and the Eqs. (2.9) and (2.10) depend on a and 
/3 only through the combination (a — 1)/ /3. Therefore, for a given w(z = 0) = wo and 
q(zj-) = 0, we can find the combination (a — l)//3 and f2 A . Note also that fi A = 1 
gives w = — 1, as expected from a universe dominated by the cosmological constant. 



It is worth examining the special case of a = 1. Integrating the Eq. (2.5) for a = 1, 



and using (2.8) to the define the corresponding EoS, gives 

2^a 



w(z) 



•1 + 



(2.11) 



3C(3 + Qn x ]og(l + z)' 

where C is the integration constant. At z = gives w = — 1 + 2Q\/(3C/3). Thus 
depending on the relative sign of /3 and C, we can have quintessence or phantom 
behavior. An interesting limit of Eq. (2.5) can be obtained when H « H 2 , which is 
equivalent to one of the slow-roll conditions used in inflationary models. In this limit 
H 2 — A/(l — a), which gives rise to exponential expansion (a ^ 1). 



3 Holographic density with matter 

In order to construct a more realistic scenario, we include the contribution of matter, 
and by using the observational data we can define the parameters of the model. 
Adding the contribution of matter, the Friedmann equation becomes 

H 2 = \ Pm0 e-^ + X + aH 2 + (3.1) 
3 2 ax 

with the solution given by (in terms of the scaled Hubble parameter) 

H 2 = Q m0 e^ + Q A + Qe- 2 ^- 1 ^ + ^"f^V 3 * (3.2) 

op — la + I 
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where Q m o = p m o/(3i?o); ^ is the integration constant and the last three terms 
represent the contribution of the holographic dark energy. The third term provides the 
possibility of crossing the phantom barrier, giving rise to quintom behavior. Note that 



the last term in (3.2) makes contribution to the matter sector, which indicates that 



the present holographic density also contributes to the presureless matter component 
of the universe, which is about 27% of the total. In what follows we will consider the 
sum of the terms proportional to e~ 3x as the matter component (in fact is an effective 



matter component) and the Eq. (3.2) will be written as 



H 1 



where 



rnCl 



3/3 - 2a + 2 



mO 



(3.3) 



(3.4) 



and Q' m0 should be adjusted with the data. According to current observations, the 



main contribution to dark energy comes from the constant term in Eq. (3.3), and the 



third term may be responsible for possible quintom behavior. The coefficients in Eq. 



(3.3) are subject to the flatness condition (at x = 0) 



n' m0 + n A + n 



(3.5) 



It is important to note here that the net contribution to the present density parameter 



of the dark energy comes from the sum of the last two terms in Eq. (3.5), which we 
can resume as 

Q DE = Q A + Q (3.6) 



so that the meaning of each separate term in (3.6) is not important, but what matters 



is their sum which should contribute about 73% to the total balance in (3.5 ). Note also 



that as in the previous case of DE dominance, a and (5 appear in the combination 
(a — !)//? in (3.3), and hence we can define the parameter 7 = (a 



l)/(3. This 

effectively reduce the number of the free parameters of the model to three, as can be 



seen from (3.3) by using the above replacement for a and (3 and the flatness condition 



(3.5). At the same time this simplifies the procedure of fitting the parameters with 
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the current observational data. In terms of the redshift z the Hubble parameter takes 
the form 

H 2 = O' m0 (l + zf + A + 0(1 + zf\ (3.7) 
and the respective total EoS becomes 



w t (z) 



-1 + 



2 7 Q(l + ^ + 3^ (l + 



(3.8) 



3 [n(i + z)^ + a^(i + z)» + i - o' m0 - o] 

where the flatness condition ( |3.5[ ) was used. Before beginning with the fixing of the 
parameters, we illustrate the behavior of the model by choosing some representative 
values of the parameters, guided by the redshift transition in the deceleration ac- 
celeration parameter and by the present observational value of the dark energy EoS 



parameter. From (3.7|) the deceleration parameter is given by 
q(z) ■■ 



-1 + 



27^(1 + ^ + 30^(1 + 



2 [0(1 + z)*i + O' m0 (l + zf + 1 - n> m0 - 0] 
The condition q(z T ) = gives O in terms of O' m0 and 7: 



(3.9) 



n 



2(1-0: 



mO) 



l + ( 7 -l)(l + z T ) 2 ^ 
Let's turn to the equation of state of the dark energy sector. 



(3.10) 



From (3.7) it follows the expression for the holographic dark energy density as 



PDE = 3H 2 (n A + n(i + z) 2 "<) 



(3.11) 



using the continuity equation pi> E + 3H(p DE + Pde) = in terms of the redshift z 
we find the DE pressure density, and the following expression for the EoS of the DE 



Pde 
Pde 



•1 + 



[1 + z) 1 dp 



DE 



3 Pde dz 



-1 + 



1 2 7 0(1 + z 



l 2 T 



3 A + 0(l + z 
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at z = (3.12) becomes 



WdEq = -1 + 



1 2 7 
31-0: 



(3.12) 



(3.13) 



mO 



where the flatness condition (3.5) was used. Assuming the current observational 
restrictions on the DE EoS [8], we can write wde = —1 + 5, where S represents the 
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deviation from —1 which may be in the range —0.2 < 5 < 0.2. Replacing in (3.13) 
we obtain the relation 
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1 2 7 1 2jQ 



3 Q A + Q 3 1 - Q' r 



(3.14) 



mO 



for a given 5 and zt, and using the relations (3.5), (3.10) and (3.14), the number of 
free parameters appearing in the solution (3.7) reduces to one. Note that for 5 = 



it follows that 7 = or Q = 0, both of them conducing to the cosmological constant. 
And negative values of 5 (indicating quintom behavior of the model), are possible if 
7 < (or even 7 > 0, < 0). 



Replacing from (3.14) in (3.10) we obtain an equation to determine {l' m0 for a given 
Zt, 8 and 7: 



36(1 - Q' m0 ) [1 + (7 - 1)(1 + zt) 2 "} = 7 [2(1 - Q' r 



mO) 



zt) 3 } (3.15) 



Thus, solving the Eq. (3.15) for 5 = 0.05, zT = 0.75 and 7 = 0.1 gives VL' m0 = 0.273. 



Quintom behavior can be obtained for 5 = —0.05 (i.e. w _d_e = —1-05), zT = 0.75 and 
7 = —0.1, giving Q' m0 = 0.269. From (3.14) Q takes the values 0.545 and 0.548 respec- 
tively, and the dark energy density parameter can be found from: Qde = + ^ = 
1 — n^Q, completing in this manner the determination of the original model parame- 



ters as displayed in (3.7). To obtain Vt' mQ we have assumed concrete values for 7, but 



we will use the observational data to fix the three parameters of the model, without 
assuming any prior values. 

Cosmological constraints from SNIa, CMB and BAO. 

We turn now to use the data sets to constrain the parameters Q' m0 , wde and 7. 
We performed a x 2 analysis to obtain the constraints on Q' m0 wde and 7 using the 
Constitution (397) SNIa data set [Sj, the shift parameter of the cosmic microwave 
background (CMB) form WMAP5 observations [UJ and the baryon acoustic oscilla- 
tions from the Sloan digital sky survey (SDSS) [TS]. The details of the x 2 technique 
and definitions of the main observational quantities used to constraint the cosmolog- 
ical parameters, can be found in [181 EJ H21 H3J HI] (see [45] for holographic dark 
energy models). We will just quote here that the luminosity distance times H for 
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our model is given by 



dr(z) 



:i+z) 



cdz' 



o H(z',n' m0 ,w DEo ,j) 



(3.16) 



where H(z, £l' m0 , W£># , 7) is given by (3.7), after replacing f2 in terms of wde through 
the relation (3.13). After minimizing the \ 2 function with respect to the model 



parameters Q' m0 , wde and 7, the following best fit values were found : f2' m0 = 0.272, 
wde = — 0.962, 7 = 0.9. This values define the best fit value for h: (the Hubble 
constant H Q in units of 100 km/s/Mpc) h = 0.649 The behavior of the Hubble 
parameter against the experimental bars (see [46] and references therein) is shown in 
Fig. 1. 



o 

Oh 




Fig. 1 The behavior of H for the best fit values Q' m0 = 0.272 7 wde — —0.962 and 
7 = 0.9, with h = 0.65. The error bars correspond to the sample of passively 

evolving galaxies 



4 discussion 

In this work we propose a new model of holographic density which incorporates a con- 
stant term as the first approximation to the infrared (IR) cut-off, and has notorious 
advantages with respect to the previous models, as it contains the ACDM model and 
additional DE term which could play an important role in the late time dynamics of 
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the universe. The cosmological constant in the context of the holographic principle 
have been discussed in [40J (where using thermodynamical arguments it was showed 
that the most probable value of the cosmological constant is zero, but a very small 
cosmological constant is still allowed), and the presence of a constant term in the 
holographic density brings a quantum gravity nature to the cosmological constant. 
Moreover, the current astrophysical data indicates that the constant cosmological 
constant adjusts very well to the observations, and this fact supports the presence 
of constant term in the holographic density. Adding this constant term to the IR 
cut-off may improve even models like the Hubble scale which was unable to produce 
accelerated expansion [32], [33]. In this case if we add a constant term to the squared 
Hubble scale, we obtain the ACDM model. In order to reduce the number of free 
parameters, we grouped the terms proportional to e~ 3x in an effective matter density 
parameter, and the dependence on the initial model parameters a and /? reduces to 
dependence on the combination 7 = (a — l)//3. By using the astrophysical data we 
can fix 7, and this means that one of the parameters a or still free. If for instance 
we let (3 free, then the relevant cosmological parameters which dictate the dynamical 
evolution of the universe do not depend on /3, and there is apparently not dynamical 
criteria to fix (5. To fix ft we can use geometrical criteria in order to convert the last 



two terms in (2.3) in a term proportional to the Ricci scalar. In this case the constant 
is automatically defined. From the relation a = 7/? + 1, if we set /3 = 1/(2 — 7) 
then the last two terms in (2.3) become 1/(2 — / y)(2H 2 + H). 

We have also shown that the holographic density contributes to the matter density 
defining in this way the "effective" matter density containing the usual baryonic and 



dark matter (see (3.2 3.4)). We can interpret this result as meaning that the present 
holographic density may account for the dark matter content, and in this sense pre- 
viously considered models of holographic DE based on local quantities H 2 and H 
can give rise only to constant DE EoS parameter. The model parameters (n^ , Q, 7) 
were constrained using the observational data including the Constitution sample of 
SNIa, the CMB shift parameter given by WMAP5, and BAO measurement from 
SDSS. We found the best fit values Q' m0 = 0.272, w DEo = -0.962 and 7 = 0.9, with 
Xmin = 465.941. The Hubble parameter versus the redshift is plotted in fig. 1. A 
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detailed study of the available data to constraint the parameters will be done else- 
where. 

Resuming, the present model could be relevant because of its simplicity; it introduces 
a "constant" cosmological term in the infrared cut-off of the holographic density, 
which puts the cosmological constant in the context of the holographic principle, at- 
tributing quantum gravity nature to the cosmological constant; it improves previously 
considered holographic models of DE based on local quantities; it may reproduce the 



standard ACDM model even without explicitly giving matter content (see (2.6) at 
(a — I)/ (3 = 3/2). We expect that future high precision experiments allow accurately 
determine the parameters of the model and define its quintessence or quintom nature. 

Acknowledgments 

This work was partially supported by Fundacion SENECA (Spain) in the frame of 
PCTRM 2007-2010. 

References 

[1] A.G. Riess, et al., Astron. J. 116, 1009 (1998), [astro-ph/9805201] 
[2] A.G. Riess, et al., Astron. J. 117,707 (1999) 
[3] S.Perlmutter et al, Nature 391, 51 (1998) 

[4] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 
517, 565 (1999) [astro-ph/9812133] 

[5] P. Astier et al, Astron. Astrophys. 447, 31 (2006) [astro-ph/05 10447]. 

[6] M. Kowalski, et. al., Astrophys. Journal, 686, p.749 (2008), arXiv:0804.4142 

[7] R. M. Amanullah et. al. , arXiv:1004.1711 [astro-ph.CO] 

[8] M. Hicken et al, Astrophys. J. 700, 1097,2009; arXiv:0901.4804v3 



12 



[9] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 
[astro-ph/0302209] 

[10] D. N. Spergel et al, Astrophys. J. Suppl. 170, 377 (2007); astro-ph/0603449. 

[11] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009), 
arXiv:0803.0547 

[12] D. Larson et al, Astrophys. J. Supp. 192, 16 (2011); arXiv:1001.4635v2 [astro- 
ph.CO] 

[13] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) [astro- 
ph/0310723] 

[14] M. Tegmark et al. [SDSS Collaboration], Astrophys. J. 606, 702 (2004), astro- 
ph/0310725 

[15] K. Abazajian et al. [SDSS Collaboration], Astron. J. 128, 502 (2004) [astro- 
ph/0403325] 

[16] K. Abazajian et al. [SDSS Collaboration], Astron. J. 129, 1755 (2005) [astro- 
ph/0410239]. 

[17] U. Seljak et al. [SDSS Collaboration], Phys. Rev. D 71, 103515 (2005), astro- 
ph/0407372. 

[18] D. J. Eisensteinet al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005) [astro- 
ph/0501171]. 

[19] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 74, 123507 (2006), astro- 
ph/0608632. 

[20] P.J.E. Peebles and B. Ratra, Rev. Mod. Phys. 75 (2003) 559 [astro-ph/0207347] 

[21] Edmund J. Copeland, M. Sami and Shinji Tsujikawa, Int. J. Mod. Phys. D 15 
1753-1936 (2006), arXiv:hep-th/0603057 



13 



[22] J. A. Frieman, M.S. Turner and D. Huterer, Ann. Rev .Astron. Astrophys. 46, 
385 (2008), arXiv:0803.0982[astro-ph] 

[23] Varun Sahni and Alexei Starobinsky, Int. J.Mod.Phys.D9, 373-444, 2000, astro- 
ph/9904398; 

[24] V. Sahni, Lect. Notes Phys. 653, 141-180 (2004), arXiv:astro-ph/0403324v3 
[25] T. Padmanabhan, Phys. Rept 380, 235 (2003), [hep-th/0212290]. 
[26] J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973). 
[27] G. 't Hooft, [gr-qc/93 10026]. 

[28] R. Bousso, JEEP 9907, 004 (1999), [hep-th/9905177]. 

[29] A. Cohen, D. Kaplan and A. Nelson, Phys. Rev. Lett. 82, 4971 (1999), [hep- 
th/9803132]. 

[30] L. Susskind, J. Math. Phys. (N. Y) 36, 6377 (1994). 

[31] W. Fischler and L. Susskind, [hep-th/9806039]. 

[32] S. D. H. Hsu, Phys. Lett. B 594, 13 (2004), [hep-th/0403052]. 

[33] M. Li, Phys. Lett. B 603, 1 (2004), [hep-th/0403127]. 

[34] E. Elizalde, S. Nojiri, S.D. Odintsov and P. Wang, Phys. Rev. D 71, 103504 
(2005), hep-th/0502082 

[35] S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38, 1285 (2006), hep-th/0506212. 

[36] L. N. Granda and A. Oliveros, Phys. Lett. B 669, 275 (2008), gr-qc/0810.3149. 

[37] L. N. Granda and A. Oliveros, Phys. Lett. B 671, 199 (2009), gr-qc/0810.3663. 

[38] L. N. Granda, Int. J.Mod.Phys.D, V.18, No.ll, 1749 (2009); arXiv:0811.4103 [gr- 
qc] 



14 



[39] L. N. Granda and A. Oliveros, arXiv:0901.0561 [hep-th] 

[40] P. Horava and D. Minic, Phys. Rev. Lett. 85, 1610 (2000), hep-th/0001145 

[41] J. R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. Roy. Astron. Soc. 291, 
L33 (1997), astro-ph/9702100 

[42] Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (2006), astro-ph/0604051 

[43] S. Nesseris and L. Perivolaropoulos, Phys.Rev.D72, 123519 (2005); arXiv:astro- 
ph/0511040 

[44] S. Nesseris and L. Perivolaropoulos, JCAP 0702, (2007) 025 ; arXiv:astro- 
ph/0612653 

[45] M. Li, X. Li, S. Wang and X. Zhang, JCAP 0906, (2009) 036; arXiv:0904.0928 
[46] D. Stern et al., JCAP 02 (2010) 008 



15 



